Chiral vibrations in the A=135 region 
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Chiral vibrations are studied with the RPA plus self-consistent tilted axis cranking formalism in 
the A=135 region. In this method chiral vibrations appear as a precursor to the static chiral regime. 
The properties of the RPA phonons are discussed and compared to experimental data. We discuss 
the limits the chiral region and the transition to the non harmonic regime. 
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I. INTRODUCTION 

Chirality is an important symmetry in many physical 
systems. In molecules it is a static property of the geome- 
try that appears in molecules with more than four differ- 
ent atoms, which can have left-handed and right-handed 
enantiomers. In these cases the left- and right-handed 
geometry is connected by reflection in a plane. In par- 
ticle physics chirality is a dynamical feature of mass-less 
particles, which indicates orientation of the intrinsic spin 
relative to the linear momentum. In a rotating nucleus 
it arises as a combination of a dynamical property, the 
angular momentum, with a static property, the reflection 
symmetric triaxial shape of the nucleus. If the angular 
momentum vector lays outside any of the three princi- 
pal planes, there are a left-handed and a right-handed 
arrangement, which are connected by the time reversal 
operation. This kind of chirality is manifested as a pair of 
(almost) degenerate rotational bands with the same par- 
ity [l| . Pairs of such bands have been seen experimentally 
in the mass regions A=105, A=135 and A=190 and the- 
oretically described using mean field tilted axis cranking 
(TAC) models @, two-particle-rotor models 0, @ or 
extensions of the IBA model [f| 0] • 

Chirality appears both in molecules and nuclei as a 
spontaneously broken symmetry: The many body Hamil- 
tonian describing these systems is invariant with respect 
to the chiral operation (space inversion or time rever- 
sal, respectively), which has the consequence that the 
exact eigenfunctions are achiral. However, there exist 
very good approximate solutions, which are chiral, where 
the left-handed and right-handed configurations have the 
same energy. The exact eigenfunctions are odd and even 
superpositions of these chiral solutions. The energy split- 
ting between these states is given by the matrix element 
between the two chiral configurations, which is the in- 
verse time for tunneling from one to the other. In most 
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molecules the tunneling time is so long that they stay in 
one of the enantiomers and the level splitting is unmea- 
surable. However, for CH 3 NHF the tunneling is rapid. 
The splitting is in the order of 100 meV, i.e. the tun- 
neling frequency is in the order of 3000 GHz, which is 
somewhat larger than the average rotational frequency 
of the molecule at room temperature. For the known nu- 
clear cases, the left-right conversion is always rapid. The 
observed energy distance between the chiral partners is 
in the order of 100 keV Q, except in very narrow spin 
regions, where the bands cross. This is comparable with 
the rotational frequency in the observed spin range. 

TAC is a microscopic mean field method that has been 
shown to very well describe the energy and the intra-band 
transition rates of the lower of the two chiral partner 
bands, see e.g. @. However, it gives either one achiral 
self-consistent solution or two degenerate chiral ones. It 
cannot describe the left-right mode, which is a well known 
deficiency of the mean field solution when it breaks spon- 
taneously a symmetry (cf. e.g. [l], U3|)- In a finite sys- 
tem, as the nucleus, the symmetry breaking develops in 
a gradual way. First a precursor appears as a slow vi- 
bration around the symmetric configuration. Becoming 
increasingly anharmonic, the vibration changes into tun- 
neling between two asymmetric configurations, which is 
progressively inhibited. 

In order to describe the splitting between the two 
bands one has to go beyond the mean field approxi- 
mation. So far this has only be done in the frame- 
work of two-particle-core coupling models [1, H, H 0]- 
These studies show a development from chiral vibrations 
toward tunneling between static left- and right-handed 
configurations with increasing angular momentum. Al- 
though the quantal nature the core-particle models ac- 
counts well for this aspect of symmetry breaking, they 
are based on several assumptions the validity of which is 
not assured. Examples are: Rigid shape and irrotational- 
flow moments of inertia for a triaxial rotor core 0, [![, 
which are not consistent with microscopic cranking cal- 
culations Q, the application of the IBA core, and the 
IB AFF coupling scheme p, 0] ■ A microscopic treatment 
starting from the TAC mean field seems important for 
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a better understanding of nuclear chirality. Obviously 
the transition from chiral vibrations to chiral rotation in- 
volves large amplitude collective motion, a treatment of 
which goes beyond the scope of this paper. However the 
regime of chiral vibration can be treated by the random 
phase approximation (RPA) approach as long as the en- 
ergy splitting between the bands is large enough for the 
assumption of a harmonic vibration remaining reason- 
able. As the energy splitting decreases the anharmonic 
effects will become more important until the transition 
point where the RPA energy becomes zero and the TAC 
mean field becomes chiral. The present paper is devoted 
to a study of the spin range before the transition point 
where the RPA approximation will allow us investigating 
the chiral vibrations and provide insight into the nature 
of nuclear chirality. 

In the static chiral picture on would expect the two 
bands to have very similar electromagnetic transition 
rates. Recent experiment have shown different intra- 
band transition rates in the two chiral partner bands in 
134 Pr @,[I3 while the bands in 128 Cs [l| and 135 Nd [H 
have similar transition rates. The different transition 
rates in 134 Pr have been interpreted as being due to cou- 
pling of shape degrees of freedom to the orientation de- 
gree of freedom p} or even that the bands are not chiral 
partner bands at all (Tlj . 

Chiral pairs of rotational bands have mainly been 
suggested for nuclei with odd proton and odd neutron 
number. In the A=135 region, they are build on the 
7r ^ l n/2 ® v ^ l i\/2 conn g ura tion. Candidate bands have 
been seen in the odd-odd nuclei 126_132 Cs, 130_134 La, 
132 " 134 Pr, 136 Pm and 138 ~ 140 Eu. In this work we fo- 
cus on the N=75 isotone chain and the Z=57 isotope 
chain which represent the central part of this region. We 
have performed calculations in the framework of tilted 
axis cranking (TAC) and random phase approximation 
(RPA) for the twin bands in the N=75 isotones and the 
Z=57 isotopes built on the {7r/i n / 2 , f/111/2} configura- 
tion. Examples for chiral partner bands exist also in 
odd-even nuclei [HI, [l4j]. The case of 135 Nd has been 
studied in [l3[ by means of the method presented in this 
paper. In principle, chiral partner bands should also ex- 
ist in even-even nuclei build on two-particle two-hole like 
configurations. For these and even more complex con- 
figurations our method can be directly applied. We de- 
scribe the formalism in section [TTJ The results for the 
energy, amplitudes and transition rates are presented in 
section IIIII An analysis of the structure of the phonons 
is given in section HVl where it will be demonstrated that 
for most of the cases the chiral character prevails, i.e. the 
coupling to the shape degrees of freedom is weak. 



II. FORMALISM 

Most earlier TAC calculations for chiral bands have 
used the Strutinsk y r enormalization for calculating the 
energies (SCTAC [15|)- Since we intend to carry out 



TABLE I: Quadrupole force strength kq. 



132 La 

134p r 

136 Pm 
138 Eu 

140 Tb 

130 La 
134 La 



kq [MeV] 



0.0367 
0.0359 
0.0348 
0.0355 
0.0363 
0.0357 
0.0372 
0.0337 



RPA calculations we need an interaction Hamiltonian. 
Hence, we choose using the self-consistent TAC with the 
QQ-force and a constant pair gap (PQTAC [15j) in two 
major N-shells (A^i ow =4 and N up =5): 



m=-2 N=4 



£ A T (P+ + P T ) - LJ ■ J 



(1) 



where ho is the spherical Woods-Saxon energy Q and 
Qm are the dimensionless quadrupole operators for each 
N-shell with an N-dependent force strength [l6| [24| : 



N low -B (2Z{N) 



N - B 



A 



1/3 



(2) 



The value of the parameter B is chosen to approxim ately 
reproduce the Strutinsky potential energy surface [la]. 
Setting the parameter B — 0.5 gives a similar prolate 
oblate mass difference and 7-softness as the SCTAC cal- 
culations in this mass region. In other mass regions a 
different value of B may be needed. The parameters 
KOi listed in table [H are adjusted to reproduce the mo- 
ment of inertia of the SCTAC calculations. The an- 
gular velocity in the cranking term, u) ■ J, is defined 
as lui — u> sin cos ip, ll>2 — w sin ?9 sin and 
U3 = ui cost?, where d and <p are the tilt angles 0. 
The monopole pair operators are denoted by P T1 and we 
use a proton and neutron pairing gap, A„ !P , that is 80% 
of the odd-even mass difference in this mass region. 

In the TAC calculations for the A=135 region, chiral 
solutions appear as a transient phenomenon where the 
angular momentum vector moves from the principal axis 
plane spanned by the 1- (short) and 3- (long) axes (ip = 0) 
to the principal plane spanned by the 2- (intermediate) 
and 3-axes (<p = 90°) [2]. The RPA method describes the 
regime when the TAC angular momentum lays in one of 
the principal planes. The excited chiral partner band can 
be described as a chiral vibration, which is a periodic mo- 
tion of the orientation of the deformed potential relative 
to angular momentum vector, i. e. a mode in the $ and 
ip degrees of freedom centered about their TAC equilib- 
rium values. There are two critical frequencies. The first 
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is associated with the transition of the TAC mean field 
solution into the chiral regime at ip = 0, and the sec- 
ond is associated with the transition out of chiral regime 
at ip = 90°. These transitions of the mean field are are 
accompanied by an instabilities of the RPA, where the 
energy of the lowest RPA phonon becomes zero. 

When the TAC solution has < \ip\ < 90° (static chiral 
regime) our present approach gives zero energy splitting 
between the two chiral configurations with ip — ±\<p\. 
However, higher order terms in our Hamiltonian would 
give rise to tunneling between the left-handed and the 
right-handed solution which causes an energy splitting 
between the two bands in this region too. Figure [T] 
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FIG. 1: Potential energy in the rotating frame for 134 Pr at 
uj — 0.4 MeV/ft. The two solid curves are calculated with 
7 or tp fixed while the dotted curve is calculated with a full 
minimization. 



shows the TAC energy in the rotating frame, spanned 
by the three principal axes of the triaxial density distri- 
bution, as a function of the tilt angle ip and the triaxial 
deformation parameter 7. The calculation is done for 
the {7r/ixi/2) vhu/2} configuration in 134 Pr at UJ = 0.4 
MeV. The energy surface is very flat especially in the <p 
direction. When we allow for the full minimization in the 
orientation degrees of freedom the surface becomes sym- 
metric in ±7, because the change of sign is equivalent 
with a reorientation of the axes. It is therefore enough 
to consider 7 > only. Figure [5] shows the potential en- 
ergy surface in the rotating frame (total routhian) as a 
function of d and (p. One can see the very soft nature of 
the potential in the ip direction. 

The RPA calculates the harmonic excitations around 
the mean field minimum. It will describe the system 
as long as we are in the chiral vibrational regime well 
before the transition to static chirality. The transition 
point corresponds to where the RPA solution goes to zero 
energy and the tilted mean field acquires a nonzero <p. 
After solving the mean field problem the Hamiltonian 
can be written as 



theta 




phi 



FIG. 2: (Color online) Potential energy surface in the rotat- 
ing frame for 134 Pr at lo — 0.4 MeV/h as a function of ■& 
and tp with constant deformation. The equipotential lines are 
separated by 25 keV. Darker color represent lower energy. 



where h m { is the diagonal mean field Hamiltonian 

h m { = E m{ + £ e kOi\a k (4) 

k 

in terms of the self-consistent mean field quasi particle 
operators a kl and H rcs is the residual interaction. We in- 
troduce the quasi-boson approximation &t = aja], where 
the 6^ are treated as exact bosons and /1 = i > j. The 
Hamiltonian H is rewritten in RPA order by only keeping 
terms up to second order in the boson operators |lOj, 

#rpa = E ml + A^by + (B^blbt + h.c.) , 



/l , V 



.(5) 

where matrices A and B are hermitian and determined 
from the residual interaction [10j]. The matrix elements 
for the QQ-interaction are 



m— — 

2 



K Ha 1v 



(6) 
(7) 



m=-2 



where g™ are the quadrupole matrix elements in quasi- 
boson, approximation (Iol | 



H — h m f + H rc 



(3) 



(8) 
(9) 



4 



and — + Bj are the two quasi-particle energies. (We 
suppressed the summation over the shell number N.) We 
solve the RPA equations 



Hkpa,0{ =E rpa 0{, 



using the strength function method of [10, Qj 
eigenmode operators Ot are 



(10) 



The RPA 



H 1 n u V 



(11) 



where the RPA amplitudes X* and are obtained by 
solving the standard set of linear equations resulting from 
Eq. (|10p together with the normalization condition 



O x ,0 



A vV 
M M 



yA^A 



#AA' 



(12) 



Since we use a separable force, this set of linear equations 
is strongly simplified [l7j |. 

There are two rotational spurious solutions in the RPA 
spectrum. One at zero energy induced by the J z operator 
and one at the rotational frequency uj induced by the J+ 
operator. Numerically the spurious solutions decouple 
from the physical RPA solutions in a stable manner if the 
mean field problem is solved accurately enough. In the 
discussion below we will only refer to the physical RPA 
solutions, which are in all cases well decoupled from the 
spurious ones. The RPA phonon energy gives the energy 
splitting between the zero-phonon lower band and the 
excited one-phonon band at a given rotational frequency 
uj. From the RPA amplitudes we derive the inter-band 
transition rates using the method of [I?) • Since RPA does 
not give any contribution to the intra-band transition 
rates, we use the TAC results [Hi for those. 



III. RESULTS 

The results for the TAC and RPA calculations for the 
N=75 isotones and the Z=57 isotopes is presented below. 
We solve the self-consistent TAC mean field problem plus 
RPA of the Hamiltonian in Eq. ([T]) for the TThn/ 2 (3vh^ 2 
configuration. 



A. Energies 

We present the results of the TAC+RPA calculations 
in Fig. [3] and [4[ which show the energy of the lowest 
RPA phonon and the experimental energy splitting be- 
tween the chiral bands, which are identified with each 
other. For reference we also plot the energy of the lowest 
two-quasiparticle excitation as well as the energy of the 
second collective RPA phonon. The nuclei studied here 
represent the chiral region for odd-odd nuclei in the N=75 
isotone chain and the Z=57 isotope chain. For lower Z or 
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FIG. 3: (Color online) Energy difference in the rotating frame 
between the two chiral partner bands in the N=75 isotone 
chain. RPA results are compared with experimental data [H, 
0, H EH, 03 ■ We also plot the second RPA phonon energy and 
the two-quasi particle energy, E2 qp . The experimental values 
show the difference between the experimental routhians of the 
two bands at grid points of to calculated by interpolation. 



larger N we approach the shell closing where the defor- 
mation disappears. Smaller N leads to disappearance of 
the triaxiality. Larger Z approaches the proton drip line, 
where very little high spin data are available. In the low- 
Z part of the A=135 region the TAC energy minimum 
sits at p = at low spin. As we move towards larger Z 
the minimum moves over to (p = 90°. With increasing lu, 
(p = 90° is preferred, because the intermediate axis has 
the largest moment of inertia 0] ■ 

In the lighter nuclides the energy of the lowest RPA 
phonon is substantially smaller than the mean field two- 
quasiparticle energy and the energy of the second RPA 
phonon. In the heavier nuclides the difference is smaller. 
For the lowest phonon being a well developed collective 
excitation it must have a substantially smaller energy 
than the next excitations. When the lowest phonon sits 
in a region where other excitations are located the collec- 
tive strength will become fragmented, and the observed 
bands will be much like two-quasiparticle excitations. 

Figure [3] shows the results for the N=75 isotones. The 
RPA results show a clear trend which reflects proper- 
ties the TAC mean field solutions. In the lighter sys- 
tems, where the TAC solutions has p = 0, the phonon 
energy decreases with spin, reaching zero at the critical 
frequency for the onset of static chirality (p > 0) . In the 
heavier systems, where the TAC solution has p> = 90° 
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already at the band head, the phonon energy increases 
with increasing angular momentum. 136 Pm is the limit- 
ing case where we have a chiral vibration around ip = 
and decreasing phonon energy at low spin, then a static 
chiral regime < ip < 90°, and finally a second chiral 
vibration around ip = 90° and increasing phonon energy. 
Comparing with experiment, 130 Cs, 134 Pr, 138 Eu have 
the right trend and magnitude of the energy splitting. 
132 La, 136 Pm have a much more constant energy splitting 
in the experiment than given by the RPA calculations. 
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FIG. 4: (Color online) Energy difference in the rotating frame 
between the two chiral partner bands in Z=57 isotopes. RPA 
results are compared with experimental data [20I . l2ll |. We 
also plot the second RPA phonon energy and the two-quasi 
particle energy, E2 qp . 

Figure d] show the results for the Z=57 isotopes. The 
magnitude of the energy splitting is reproduced but there 
is a problem with the angular frequency trend in 132 La. 
In experiment 132 La looks similar to 130 La while in the 
calculations it looks more like 134 La. 

Our RPA calculations give the correct energy splitting 
away from the region where the RPA energy approaches 
zero, which is where we can expect the RPA to work well. 
However, spin dependence is not always well reproduced. 
The RPA phonon energy tends to approach zero faster 
and more often than in experiment, which indicates that 
anharmonic effects are needed to understand the data 
in this regime. The RPA phonon energy is generally in- 
creasing or decreasing with rotational frequency. The ex- 
perimental energy splitting, on the other hand, does not 
change much with u> in some cases, which seem to ap- 
pear predominantly in the frequency regions where the 
TAC gives static chirality, < ip < 90°, i.e. TAC+RPA 
predicts zero splitting. The constant energy splitting in 
these regions should be attributed to the tunneling be- 
tween the left- and right-handed solutions. Its relatively 
large value of about 300keV indicates strong mixing be- 
tween the two chiral configurations, which is consistent 



with the flat potential in <p direction in Fig. [2] and our 
analysis of the composition of the RPA phonon wave 
function discussed in section HVl below. In 134 Pr the two 
band cross over in experiment. In 130 Cs and 135 Nd (not 
shown here, see [IH) one sees in experiment an avoided 
crossing with some interaction. In all three cases, the 
crossing frequency correlates with the frequency where 
the phonon energy goes to zero. At the moment, it re- 
mains unclear why these nuclei are different from 132 La 
and 136 Pr with a nearly frequency-independent energy 
difference between the bands. To understand this phe- 
nomenon further it is clearly necessary to take into ac- 
count large amplitude anharmonic effects, which will be 
done in a forthcoming publication (22j. 



B. Transition rates 

In addition to the energy splitting, the RPA also gives 
the transition rates between the two chiral partner bands. 
In Fig. rSlTUl we show the calculated B(E2) and B-(Ml) 
inter-band rates. The intra-band values calculated by 
means of TAC 
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are also included. Within the present 
approach they are the same for the two bands. The inter- 
band reduced transition probabilities are given by 



i) = (l\e T r 2 Y^\o) 



-,2 
(13) 



B(E2,I^I 
and 

B(Ml,I->I-m) = (l l/i 1 ^ 0) 2 m = -1,0,1, (14) 

where |0) is the RPA ground state and |1) is the first 
excited RPA state. The charges e T are for 1 protons and 
neutrons. The spin parts of the magnetic operator are 
attenuated by the factor of 0.7. We include terms up to 
linear RPA order in the quadrupole, r 2 Y"2m and magnetic 
fj, m transition operators. The transition operators are 
transformed from the intrinsic frame into the lab frame 
by means of the relation 



rplah 
1 lm 



(15) 



where the —n/2 rotation brings the intrinsic system from 
x-axis to the z-axis representation. 

The inter-band Ml transitions are much stronger in the 
lighter N=75 isotopes where the B(M1) values are typi- 
cally around 1/3 of the intra-band values. In the heavier 
N=75 isotopes as well as in 134 La the inter-band B(M1) 
are weak. For the E2 transitions in the lighter N=75 
isotopes the dominant inter-band transition is I — » 1—2, 
except approaching the critical frequency the transition 
I — ► I ± 1 increases in strength. In the heavier N=75 
isotopes B(E2,I — ► I + 1) becomes the largest inter- 
band value. However, the J^/+lor/ + 2 will often 
not be seen in experiment due to their low or negative 
energy. The ratios of the different transition proba- 
bilities vary substantially over the studied region. The 
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FIG. 5: The RPA inter-band transition rates in 130 Cs, 132 La 
and 134 Pr. 



FIG. 7: The RPA inter-band transition rates in 136 Pm, 138 Eu 
and 140 Tb. 
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FIG. 6: The TAC intra-band transition rates in 130 Cs, 132 La 
and 134 Pr. 
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FIG. 8: The TAC intra-band transition rates in 136 Pm, 138 Eu 
and 140 Tb. 



B{E2 1 I -> I± 1) and B(M1,I -> I± 1) are typically 
the strongest. 

There are experimental data for the inter-band transi- 
tion rates in 134 Pr [23]. Our calculated rates are some- 
what larger than what is seen in experiment. There are 
also data available for the ratio of B(Ml out )/ B(Ml in ) in 
i36p m |JJ| ^ unfortunately not in the frequency range 
where we have RPA results. 



IV. SHAPE AND ORIENTATION AMPLITUDES 

The energy of the RPA phonon and the transition rates 
do not give direct information about the degree the RPA 
phonon is dominated by orientation or shape fluctua- 
tions. The potential energy surfaces as Fig. [T] and [2] give 
some hint of which are the important degrees of freedom 
but they do not provide conclusive information, because 
the missing collective mass parameter is as important as 
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FIG. 9: The RPA inter-band transition rates in 130 La and 
134 La. 
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the curvature of the potential. In order to analyze the 
nature of the RPA phonons we calculated the fluctua- 
tions of the five quadrupole operators, which are then de- 
composed into fluctuations of the Euler angles and shape 
fluctuations. 

We define the fluctuation of the quadrupole operator 
Qm generated by the RPA phonon as the difference of the 
dispersions in the one-phonon state |1) = 0[ |0) and zero- 
phonon state |0), where 0\ creates the lowest phonon, i.e. 



<l|(Q ro -(l|Q ro |l)) 2 |l>-(0|(Q r , 



(16) 

<0|QJ0» 2 |0). 



Using that in RPA approximation (1 \Q m \ 1) = 
(0 \Q m \ 0) = q m wc find 



AQ 2 m = 2[Q m ,0\ 
We define, AQ m , for m — —2, 2 as 
AQ m = V2 \Q m ,0\ 



(17) 



(18) 



where the over all sign is ambiguous due to the freedom of 
the phase of the wave-function, but the relative signs be- 
tween the different components are relevant. The body- 
fixed frame is attached to the TAC solution, the shape 
of which of is given by the two self-consistent mean field 
quadrupole moments go — {Qo) and q2 = q-2 = {Qv)- 
Its orientation in the laboratory frame is given by the 
three Euler angles ip, ip. We decompose the AQ m 
into a part that is generated by a reorientation of the 
body- fixed frame with constant q m , which is expressed 
in terms of changes of ip, i?, ip, and into a part that is 
generated by changes of the shape, which are expressed 
in terms of changes of the deformation parameters e and 
7. By means of the Wigner functions D^ nm ,{ijj,'d,ip), 
we express the five components of the quadrupole op- 
erators Q m in the laboratory system in terms of the 
two intrinsic quadrupole moments Q = Qe cos 7 and 
Qi = Q-2 = Qesin7/%/2, i.e. 

Q m (iP, V, if, e, 7) = D 2 m0 Q + (D 2 m2 + D 2 m _ 2 ) Q 2 , (19) 

where Q is a constant, e, 7 are the standard deforma- 
tion parameters, and the three Euler angles, tp, and 
ip specify the orientation of the body-fixed coordinates 
in the laboratory frame. We derive linear expressions 
for AQ m in terms of the fluctuation in Euler angles and 
shape parameters 7 and e 



dot; 



eiM^e}. (20) 



Inverting this system of linear equation, we express the 
fluctuation of Euler angles and of the shape parameters 
in terms of AQ rn and mean field expectation values, 



Aip = 
Atf = 

Aip = 

Ae 

e 

A7 = 



AQi 

2q 2 ' 
AQ-i 

q 2 + V2qo ' 

Ag_ 2 

q 2 - \/2q ' 

(V3g 2 - g0 ) AQp + (V3g 2 + go) AQ 2 

2 (q 2 + q\) 
(\/3g 2 + go) AQo - (\/3g 2 - g ) AQ 2 

2(<Zo 2 +9 2 2 ) 



(21) 
(22) 
(23) 

, (24) 
• (25) 



The results are shown in Fig. [TT] and [T2J One can 
clearly see that the lowest phonon in this region is dom- 
inated by orientation fluctuations especially in the ip di- 
rection, which is consistent with the TAC potential en- 
ergy surface in Fig. [2] The amplitude of Aip is typically 
30-40°. The fluctuations in tp diverge at the point where 
the RPA energy goes to zero, i.e. at the transition to 
aplanar chiral rotation. The fluctuations in d are 10- 
20° while the fluctuations in ip are half that. In all these 
cases, the shape fluctuations are small in amplitude com- 
pared to the orientation fluctuations. The fluctuations 
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FIG. 11: (Color online) The fluctuation of Euler angles in the 
RPA phonons in the N=75 isotone chain. Cs, La and Pr in 
the left panels and Pm, Eu and Tb in the right panels. 



70 
60 
50 
9- 40 

< 30 
20 
10 


60 

40 

CD 

< 20 



-20 
60 

40 

< 20 



- 20 o. 



- 1 


1 


1 - 
/ 

S 
s _ 






— h 1 ^ 


— 1 — 
'-^4 — 


— i — 

130, 

La 
" 2 La 

l"La._,- 






1 

i 


h^T- 

1 



0.2 0.3 
hco [MeV] 



0.4 



FIG. 13: (Color online) The fluctuation of Euler angles in the 
RPA phonons in the Z=57 isotones. 



- 1 1 \ 

" / ™Cs" 

/ ""La _ 

' 4 Pr - 

1 1 -1 1 1 


I 1 1 1 1 1 1 

156 Pm- 

l58 Eu " 

""Tb" 

-v - 

— 1 1 1 1— 


- 1 S. 1 

- / • - 

/ 

/ 

-,1,1,1,- 


" , 1 , 1 , " "l"^^ 



0.1 



0.2 0.3 0.4 
hco [MeV] 



0.1 



0.2 0.3 0.4 
hco [MeV] 



0.5 



5 



80 
60 
40 
20 

-20 
-40 
-60 

10 

5 




-1Q. 



- 1 ; 
~' * . \ 


i 


" 2 La 

134 

■■■ La .l ~ 


I 





0.1 



0.2 0.3 
hco [MeV] 



0.4 



FIG. 12: (Color online) The fluctuation of shape parameters 
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in e are few percent of the equilibrium values and the 
fluctuations A7 < 6° for all cases but 134 La (cf. next 
paragraph). This is in contrast to Ref. 0, which, taking 
into account the shape degrees of freedom in the frame- 
work of the the IBAFF two-particle-core approach, found 
a strong coupling between the shape and orientation de- 
grees of freedom. In our microscopic approach, the lowest 
phonon is a rather clean chiral vibration. 

Figures [T3] and [TJ] show how the structure of the low- 
est RPA phonon develops when moving away from the 
N=75 chain. In 130 La it is dominated by orientation 
fluctuations but not in the tp direction as for the N=75 
cases. This nuclide also has relatively strong 7 fluctua- 
tions especially at large u>. In 134 La with N=77 we are 
approaching the N=82 shell closure. This leads to a re- 
duction in the deformation and increased amplitudes of 
the shape fluctuations in the RPA phonon. The case of 



FIG. 14: (Color online) The fluctuation of shape parameters 
in the RPA phonons in the Z=57 isotones. 

134 La can no longer be thought of as a chiral vibration 
even though it has large components of orientation fluc- 
tuations in the wave function. The large jump in several 
of the 134 La curves is due to the level crossing between 
the first and second RPA solution seen in Fig. 0J 



V. CONCLUSIONS 

We have used the RPA method to study the energy 
splitting between chiral partner bands. The RPA phonon 
energy gives the correct size of the energy splitting at low 
spin in most cases but the trend with increasing angular 
momentum is not always reproduced. The RPA calcula- 
tions also give the inter-band transition rates. As the en- 
ergy splitting between the chiral partner bands becomes 
small the anharmonic effects are expected to get more 
and more important. 
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By analyzing the RPA amplitudes we found that near 
Z=57 and N=75 the lowest phonon is a collective chiral 
vibration. It is dominated by large amplitude fluctua- 
tions in the orientation angle (p. The amplitudes of the 
shape oscillations are found to be much smaller than the 
ones of the orientation oscillations. The large amplitude 
nature of the (^-fluctuations indicate that for a full under- 
standing of the chiral bands one has to take into account 
effects that go beyond the harmonic RPA approximation 
in the orientation degrees of freedom. Anharmonicitics 
become dominant where the energy of the first RPA so- 
lution goes to zero and the fluctuation amplitude in (p 
reaches 90° . Studies of the consequences of anharmonic- 
ities for the nature of the band crossing and the source 
of the different in band transition rates seen in 134 Pr are 
to be published in the near future [22| . 

The quite pure chiral nature of the RPA solution 
seems to be localized around the N=75 isotones. Mov- 



ing to more neutron-rich nuclei structure of the first RPA 
phonon changes to a complex mixture of shape and orien- 
tation degrees of freedom, reflecting the approach of the 
shell closure at N=82. In the more neutron deficient nu- 
clei the e deformation is increasing and the 7 deformation 
is decreasing which reduces the collectivity of the lowest 
RPA solutions that takes on two-quasiparticle character. 

The calculation of the inter-band transition rates show 
that the ratios of the different Al components of the 
transitions vary relatively fast with N and Z. 
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